AD-A149  336  SURFACE  STRUCTURE  ANALYSIS  USING  LOH  ENERGY  ION 
SCATTERING  I  CLEAN  FECOOl.  .  (U)  PENNSYLVANIA  UNIV 
PHILADELPHIA  DEPT  OF  MATERIALS  SCIENCE  AND  E.  . 
UNCLASSIFIED  L  MARCHUT  ET  AL.  1984  TR-28  F/G  7/4 


jOPY  AD-A149  336 


REPORT  DCC'JMEHTATIOM  PAGE 


►  KAI)  IV.THC'  Tr'i*.-. 
nri-riM-  I~f  w<r\TV.r.  |  r  if 


UOVT  ACCESSION  «0.|  J  BEClRltNT's  CAT 


Technical  Report  No.  28 


«  TlTct  f An*  J.6IITI.) 


Surface  Structure  analysis  using  low  energy  ion 
scattering  I.  Clean  Fe(001) 


*  PyfcRlOO  COvCftCO 


(.  PCR*o*MiMG  one.  report  numicr 


7.  AuTmOM'ij 

L.  Marchut,  T.M,  Buck,  G.H.  Wheatley  and  C.J. 
McMahon,  Jr.  : 


i.  contract  or  ciunt  NUMatx'n 


N00014-83-K-0565 


*  rerforminS  organization  name  and  aooresj 

Trustees  of  the  Univ.  of  Pennsylvania 
Office  of  Projects  Research  and  Grants 
3451  Walnut  St. 

Phi  1  acIpI  nhia 


>1  CONTROLLING  Office  NAME  anO  ADORESS 

Office  of  Naval  Research 
Dept,  of  (he  Navy 
500  N.  Quincy  St. 


10.  BROG*am  ElEmEnt  »»SjECT.  Ti'i« 
ABE  A  A  WORK  unit  NUMBERS 


NR  629-844/6-08-84  (410) 


11.  RERORT  CATE 

November  1984 


IS.  NUMBER  Of  PAGES 


u  -  agency  name  *  ADDRESS/;/  dillatani  from  Controlling  Oltic a)  IS.  SECURITY  CLASS,  fof  thf  tapertj 

Unclassified 


1 1  w.  declassification  oownoRa;  ns 
SCHEDULE 


|1S.  OiSTRiSuTiOn  STATEMENT  'at  rhi,  fi.porij 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (ot  thb  •  detract  antarad  in  Block  20,  ll  dilfarant  from  Report) 


JAM  1  0  1985 


it.  supplementary  notes 


V  A 


IS.  KEY  WORDS  (Continue  on  rivana  aid*  H  nacaaia/y  and  t  dart  (Ur  by  block  numbar) 


20.  ABSTRACT  (Continue  on  revert#  aid a  If  nacaaamrr  and  Identify  by  block  numoar) 

^  A  clean  Fe(001)  unreconstructed  surface  has  been  analyzed  using  low  energy 
ion  scattering  with  time  of  flight  energy  analysis.  The  intensity  of 
scattered  9.5  KeV  Ne+  ions  and  neutrals  was  measured  as  a  function  of 
elevation  and  azimuth  of  the  incoming  beam  while  keeping  a  fixed  scattering 
angle  of  90°.  The  observed  intensity  variations  can  be  described  by 
considering  the  shadowing  and  blocking  of  the  second  and  third  atomic  layers. 

A  simple  Firsov  potential  was  used  to  calculate  the  shadow  cone  dimensions.  The 
good  agreement  between  the  observed  and  calculated  variatiuns  tn -scattered  _ 


DD  1473  COITION  OF  1  NOV  it  It  OBSOLETE 

S/N  OlOJ-OM-AtOI  i 


Security  cl aj»ific ation  of  thii  race  i>.i.  *«*.» •<»> 


84  12  31  -  ob'r 


•elativel 


simple  technique 

\ 


UNCLASSIFIED 


Surface  Science  141  (1984)  $49-566 
North-Holland,  Amsterdam 


549 


SURFACE  STRUCTURE  ANALYSIS  USING  LOW  ENERGY  ION 

SCATTERING 

I.  Clean  Fe(001) 

L.  MARCHUT,  T.M.  BUCK,  G.H.  WHEATLEY  and  C.J.  McMAHON,  Jr.  * 

AT&T  Bell  Laboratories,  Murray  Hill,  New  Jersey  07974,  USA 
Received  9  August  1983;  accepted  for  publication  28  February  1984 


A  clean  Fe  (001)  unreconstructed  surface  has  been  analyzed  using  Low  Energy  Ion  Scattering 
with  lime  of  flight  energy  analysis.  The  intensity  of  scattered  9.5  KeV  Ne  *  ions  and  neutrals  was 
measured  as  a  function  of  elevation  and  azimuth  of  the  incoming  beam  while  keeping  a  fixed 
scattering  angle  of  90°  The  observed  intensity  variations  can  be  described  by  considering  the 
shadowing  and  blocking  of  the  second  and  third  atomic  layers.  A  simple  Firsov  potential  was  used 
to  calculate  the  shadow  cone  dimensions.  The  good  agreement  between  the  observed  and 
calculated  variations  in  scattered  intensity  demonstrates  the  utility  of  this  relatively  simple 
technique  for  preliminary  analysis  of  surface  structure. 


1.  Introduction 

The  Low  Energy  Ion  Scattering  (LEIS)  technique  utilizes  a  beam  of  mono- 
energetic  ions  as  a  probe  of  the  surface  and  can  be  used  for  elemental 
identification,  quantitative  analysis  and  surface  structure  analysis.  When  a 
beam  of  accelerated  ions  strikes  a  crystal  surface,  many  of  the  ions  make 
elastic  collisions  with  target  atoms.  The  energy  of  an  ion  making  such  an 
elastic  collision,  here  expressed  as  the  ratio  of  final  energy.  £,  to  initial  energy. 
£0  can  be  derived  using  energy  and  momentum  conservation  laws  as: 

£,/£u  -  [  M* /( M,  +  M,  )2]  [cos  «L  +  (  Ml /Ml  -  sin2«L )'  J]\  (1) 

Here  A/,  is  the  mass  of  incoming  ion.  M2  the  mass  of  the  target  atom  and  0L 
the  laboratory  scattering  angle.  By  measuring  the  energy  distribution  of 
scattered  particles  the  elements  present  in  the  region  sampled  by  the  ion  beam 
can  be  determined  in  a  straighforward  manner.  This  is  illustrated  in  fig.  1 
which  is  a  spectrum  for  9.S  keV  Ne  *  ions  and  neutrals  scattered  from  a  Fe 
surface  containing  segregated  Sn.  (The  Fe-3.9at%Sn  single  crystal  used  for  this 
study  could  be  sputtered  and  annealed  at  500-600  °C  resulting  in  an  ordered 
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Fig.  1.  Typical  Time-of-Flighi  LEIS  energy  spectrum  derived  from  TOF  distnbution  9.5  kev  Ne 
scattered  from  an  Fe-Sn  single  crystal.  (001)  surface,  with  segregated  Sn. 


surface  layer  enriched  in  Sn,  or  sputtered  and  annealed  at  400  °C  to  remove 
damage  due  to  sputtering  but  allow  no  Sn  segregation.  The  analysis  performed 
on  the  surface  enriched  in  Sn  will  be  the  subject  of  a  future  paper.)  Peaks  are 
seen  at  E/E0  -  0.43  and  0.70  due  to  scattering  through  90  °  from  Fe  and  Sn. 
This  energy  spectrum  was  derived  from  a  time-of-flight  (TOF)  spectrum  which 
was  calibrated  with  reference  to  90°  scattering  of  Ne*  from  Sn  using  a 
Tennelec  TC850  time  calibrator.  The  discrepancies  in  the  observed  peak 
positions  and  those  calculated  from  eq.  (1)  probably  arise  from  inelastic  energy 
loss  and  systematic  errors  in  the  time  calibration.  Such  discrepancies  are 
unimportant  to  the  analysis  presented  here.  The  background  in  this  figure  is 
from  ions  which  scattered  from  more  than  one  target  atom  before  leaving  the 
crystal. 

The  LEIS  (TOF)  technique  used  in  this  study  allows  detection  of  both  ions 
and  neutrals,  avoiding  complications  due  to  neutralization  effects  [1.2]  which 
must  be  considered  if  an  electrostatic  analyzer  (ESA)  is  used.  A  weakness  of 
the  LEIS  (ESA)  technique  is  that  many  ion-target  collisions  will  result  in 
displacement  and  removal  of  surface  atoms.  The  ion  dose  required  by  the  LEIS 
(TOF)  technique  is  100-1000  times  smaller  than  for  the  LEIS  (ESA)  method 
using  noble  gas  ions  for  two  reasons:  (a)  both  scattered  ions  and  neutrals  are 
collected  and  (b)  LEIS  (TOF)  is  a  multi  channel  technique  collecting  at  all 
energies  concurrently  rather  than  sequentially  as  with  an  ESA.  At  low  energies 
the  scattering  cross  sections  are  relatively  large,  allowing  more  rapid  data 
accumulation  than  is  possible  at  high  energies.  In  addition  the  large  scattering 
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cross  sections  lead  to  strong  attenuation  of  the  incoming  and  outgoing  particles 
so  that  analysis  can  be  limited  to  the  top  atomic  layer  [3]. 

In  this  and  many  previous  LEIS  studies  [4-11]  information  about  the 
arrangement  of  adsorbed  atoms  on  the  surface  has  been  obtained  by  taking 
advantage  of  the  shadowing  effect.  This  occurs  because  ions  approaching  a 
target  atom  are  deflected  away  from  it,  with  the  ions  approaching  more  closely 
being  deflected  the  most.  As  a  result  the  beam  does  not  penetrate  a  region 
called  the  shadow  cone  behind  each  target  atom  as  shown  in  fig.  2.  By 
comparing  changes  in  measured  scattered  intensity  with  the  projections  of 
shadow  cones  on  different  atoms,  the  arrangement  of  atoms  on  the  surface  can 
be  deduced. 


2.  Experimental  procedure 

The  LEIS  measurements  were  performed  in  a  system  described  previously 
[1,12-14],  a  diagram  of  which  is  shown  in  fig.  3.  The  cryo-pumped  UHV 
chamber  reached  an  ultimate  pressure  of  about  2x10' 10  Torr.  with  the 
pressure  increasing  to  about  2  x  10 Torr  due  primarily  to  Ne  during  the  ion 
scattering  measurements.  The  chamber  was  equipped  with  a  LEED-AES 
system  and  an  ion  gun  (both  from  Physical  Electronic  Industries.  Inc.),  in 
addition  to  the  ion-scattering  facilities. 

The  sample  was  mounted  on  a  (Varian)  heater  and  sample  temperature  was 
measured  using  a  chromel-alumel  thermocouple  spot-welded  to  a  Ta  clamp 
which  held  the  sample  onto  the  heater.  The  heater  and  sample  were  mounted 
on  a  manipulator  that  allowed  three  translational  motions  with  a  reproducibil¬ 
ity  of  ±  1  mm  and  rotation  about  two  axes  accurate  to  ±  1  ° . 


r 


Fig.  2.  Sketch  showing  the  general  relationship  between  the  trajectory  of  scattered  particles  and 
impact  parameter  p,  and  the  formation  of  a  shadow  cone. 
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Fig.  3.  Diagram  of  (he  system  used  for  LEIS  measurements. 


The  ion  scattering  was  performed  using  a  pulsed  9.5  keV  beam  of  Ne  +  ions. 
Scattered  particles,  both  ions  and  neutrals,  were  collected  by  a  Johnston 
Multiplier  -  1.0  m  from  the  sample  and  oriented  at  a  scattering  angle  6  of 
90  ° .  Time-of-flight  energy  analysis  of  scattered  particles  was  used  to  that  both 
ions  and  neutrals  were  collected,  avoiding  complications  due  to  neutralization 
effects  mentioned  earlier.  The  time  resolution  of  the  system  was  about  60  ns. 
For  9.5  keV  Ne*  scattering  from  Fe  at  4.47  keV  and  4.5  fis  time  of  flight  this 
corresponds  to  an  energy  resolution  of  0.3  keV  or  1E/E  -  0.026. 

The  sample  studied  was  from  an  Fe-1.65at%Sn  alloy  which  had  been 
directionally  solidified,  resulting  in  large  collumnar  grains  which  were  homoge¬ 
nized  at  950 °C  for  75  h  *.  A  single  crystal  with  a  (001)  surface  was  spark-cut 
from  the  ingot.  An  optically  flat  surface  was  prepared  by  mechanical  polishing 
using  alumina  powders,  followed  by  Syton  (TM)  polishing  (0.01  |im  alumina 
suspended  in  NaOH)  and  etching  with  dilute  nitric  and  glacial  acetic  acid 
leaving  a  damage  free  surface. 

The  crystal  surface  was  cleaned  further  inside  the  vacuum  chamber  by  argon 


Alloy  prepared  through  the  courtesy  of  Dr.  C.L.  Brunt  of  the  General  Electric  Research  and 
Development  Center  Schenectady. 
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ion  sputtering  (0.5  kV,  10  ^A)  with  the  sample  at  400°C.  followed  by  a  30  min 
anneal  at  400  °C  to  eliminate  damage  due  to  sputtering.  AES  analysis  of  the 
sample  at  this  point  revealed  the  presence  of  some  C  and  N  which  apparently 
had  segregated  to  the  surface  during  the  400  °C  anneal  and  could  not  be 
completely  removed  by  successive  sputtering  and  annealing.  The  LEED  pat¬ 
tern  observed  was  essentially  p(l  x  1)  with  faint  split  1/2.  1/2  spots,  probably 
due  to  C  and  N  on  the  surface.  A  small  amount  of  Sn  also  observed  in  both 
AES  and  LEIS  spectra  after  the  treatment  described  above.  The  Sn  concentra¬ 
tion  was  estimated  to  be  3.9  at%.  representing  only  a  slight  enrichment  of  the 
surface. 


3.  Results  and  discussion 

After  cleaning  of  the  sample,  LEIS  spectra  were  recorded  for  different 
orientations  of  the  beam  relative  to  the  sample.  The  coordinates  used  to 
describe  the  beam  orientation  are  shown  in  fig.  4.  The  elevation  angle  of  the 
beam  relative  to  the  sample  surface  is  designated  with  ^  =  90  normal  to  the 
surface.  The  azimuthal  angle,  <f>.  is  the  angle  of  rotation  about  the  surface 
normal,  with  </>  =  0  being  parallel  to  a  (110)  plane.  For  clarity,  when  the  beam 
is  aligned  parallel  to  a  low  index  plane  of  the  crystal,  the  indices  of  this  plane 
will  be  used  to  define  the  azimuth. 

If  relative  intensities  were  compared  using  peak  heights  without  background 
subtraction,  the  change  in  background  due  to  multiple  scattering  (which  is 
large  for  some  orientations)  would  give  spurious  results.  In  addition  the  peak 
height  may  not  accurately  represent  single  scattering  from  deeper  layers  which 
primarily  broadens  the  peak.  To  avoid  these  problems,  relative  intensities  were 
compared  using  integrated  peaks  after  subtracting  the  multiple  scattering 


Fig.  4.  Coordinates  used  to  describe  scattering  of  an  ion  by  a  surface  atom,  y  is  the  elevation  of  the 
incoming  beam  above  the  plane  of  the  surface.  *  is  azimuthal  angle  or  angle  of  rotation  about  the 
surface  normal  with  respect  to  some  reference  plane  and  0  is  the  scattering  angle. 
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background.  The  area  integrated  was  that  within  the  single  scattering  peaks, 
with  the  subtracted  background  being  the  average  of  the  number  of  counts  at 
the  two  limits  of  integration,  multiplied  by  the  width  of  the  peak.  The  single 
scattering  peaks  shown  in  fig.  1  (which  represents  a  (110),  35°  beam  orienta¬ 
tion)  have  clear  limits;  the  transition  between  the  peak  and  background  is 
sharp  and  a  small  change  in  the  limits  of  integration  would  not  effect  the  sum 
appreciably.  However,  for  some  beam  orientations  this  was  found  not  to  be 
true.  Fig.  5  shows  a  spectrum  from  clean  Fe  with  the  beam  at  (210).  45° 
(azimuth  parallel  to  (210)  planes,  elevation  45  0 ).  The  multiple  scattering  yield 
is  very  large,  and  the  limits  of  the  single  scattering  peak  are  not  obvious.  The 
limits  of  integration  for  peaks  such  as  this  have  been  chosen  to  be  reasonably 
close  to  those  used  for  the  sharp  single  scattering  peaks,  but  it  is  obvious  that 
this  integration  will  have  some  inherent  error.  The  broad  shoulders  on  the  low 
and  high  energy  sides  of  the  Fe  single  scattering  peak  in  figs.  1  and  5  are 
probably  due  to  multiple  scattering  sequences  involving  several  atomic  layers 
which  neutralize  the  ions  quite  completely  so  that  such  structure  does  not  show 
up  very  strongly  in  ESA  spectra.  Identification  of  such  sequences  has  been 
made  by  computer  simulation  for  another  system.  Cu}Au(100)  (15,16]. 

A  correction  must  be  made  to  the  integrated  intensities  measured  at 
different  elevations.  The  spot  size  of  the  beam  on  the  sample  depends  on  the 
elevation  angle;  the  intensities  were  multiplied  by  sin(^)  to  correct  for  this. 

Sequences  of  spectra  were  recorded  at  -  4°  intervals  of  rotation  of  the 


0.5  1.0 


E/E. 

Fig.  5.  LEIS  spectrum  from  the  clean  surface  with  the  beam  oriented  parallel  to  (210)  planes. 
^-45°. 
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crystal  about  its  normal  at  fixed  elevation  angles,  or  at  2.5-5°  steps  of 
elevation  with  fixed  azimuthal  angles.  The  data  from  the  azimuthal  and  tilt 
scans  have  been  combined  in  plots  of  scattered  intensity  versus  ip  and  <t>,  which 
H  will  be  called  intensity  plots  (e.g.  fig.  6). 

Ideally,  these  data  would  be  plotted  as  contour  maps  with  contours  of 
different  scattered  intensity  shown,  but  measurements  were  not  made  at 
enough  orientations  to  do  this  accurately.  The  intensity  plots  show  the  sum  of 
each  single  scattering  peak  multiplied  by  sin(  <p )  and  rounded  off  to  the  nearest 
thousand  counts.  This  averages  out  fluctuations  due  to  experimental  error  and 
JJJ  allows  the  matching  up  of  intentities  where  the  scans  intersect.  Thus  the 

intensity  plots  show  only  the  larger  intensity  variations  and  not  the  finer 
details  of  the  same  order  of  magnitude  as  the  experimental  error. 

The  intensity  distributions  were  plotted  for  *p  from  0°  to  90°  and  from 
<t>  =  0°  (a  (110)  azimuth)  to  <f>  =  45  °  (a  (100)  azimuth).  The  limits  reflect  both 
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the  possible  range  of  ^  for  90°  scattering  and  the  range  of  </>  necessary  to 
represent  the  bcc  (001)  surface.  The  intensity  plot  for  the  clean  Fe  surface  is 
shown  in  fig.  6  and  comprises  data  from  azimuthal  scans  at  \p  =  20  ° ,  35  0  and 
45°  and  from  tilt  scans  along  (110),  (310)  and  (100)  planes.  The  lowest 
intensities  ( -  1000  counts)  are  seen  near  the  (110)  azimuth  where  only  the  first 
layer  is  visible  to  the  beam.  Apparently,  single  scattering  from  the  second  layer 
does  not  occur  at  any  elevation  along  this  azimuth.  This  is  in  agreement  with 
the  shadow  cone  analysis  to  be  discussed  below.  The  single  scattering  intensi¬ 
ties  are  roughly  twice  as  large  near  the  (100)  azimuth,  showing  little  change 
with  elevation.  The  highest  intensities  are  seen  near  the  (310)  azimuth  at 
elevations  between  33  and  60  degrees. 

This  intensity  plot  will  be  discussed  in  the  following  section  where  the 
changes  in  scattered  intensity  are  compared  with  changes  in  the  number  of 
atoms  visible  to  the  beam  as  predicted  by  a  shadow  cone  analysis. 


4.  Shadow  cone  analysis 

As  mentioned  previously,  when  a  beam  of  ions  is  incident  on  an  atom,  there 
will  be  a  region  beyond  the  atom  into  which  the  ions  cannot  penetrate.  The 
formation  of  this  region  can  be  understood  with  the  help  of  fig.  2.  A  target 
atom  is  shown,  and  trajectories  of  scattered  ions  are  shown  schematically.  The 
scattering  angle  and  energy  loss  of  an  ion  due  to  a  scattering  event  can  be 
expressed  in  terms  of  an  impact  parameter,  p,  which  is  the  distance  of  closest 
approach  the  ion  and  target  atom  would  have  if  no  scattering  occurred.  The 
ions  with  small  impact  parameters  are  scattered  through  large  angles  while  ions 
with  large  impact  parameters  are  deflected  only  slightly.  The  net  result  is  that 
the  ions  do  not  penetrate  a  region  beyond  the  target  atom.  If  the  scattering 
angle  is  known  as  a  function  of  impact  parameter,  the  dimensions  of  this 
region,  called  the  shadow  cone,  can  be  calculated. 

The  shadow  cone  dimensions  used  here  have  been  calculated  in  the  manner 
of  Martynenko  [17,18]  who  applied  a  shadow  cone  analysis  to  single  crystal 
sputtering.  The  calculations  are  fairly  simple  for  the  interaction  potential  used, 
which  is  of  the  form  V{x)  -  A/x2,  where  A  includes  the  part  of  the  screening 
function  which  is  independent  of  x.  If  the  potential  is  used,  and  it  is  further 
assumed  that  near  the  edge  of  the  shadow  cone  the  ions  are  scattered  through 
small  angles  since  the  ions  and  target  atom  do  not  approach  closely,  the 
scattering  angle.  0.  and  impact  parameter,  p.  are  related  by  a  simple  analytical 
expression.  Using  this  relation,  Martynenko  derived  the  following  formula  for 
the  radius,  r,  of  the  shadow  cone  at  a  distance.  /,  beyond  the  target  atom. 

r=2.2(IA/(p+  1)£0)'/5. 


(2) 
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Here  ju  is  the  ratio  of  the  mass  of  the  incident  ion  to  that  of  the  target  atom. 

The  potential  used  by  Martynenko  is  the  approximation  to  the 
Thomas-Fermi  potential  by  Firsov  [191.  Firsov  showed  that  in  the  energy 
range  1-50  keV  the  following  screened  coulomb  potential  is  within  10%  of  the 
Thomas-Fermi  potential: 


U(r) 


Z\  Z^e 2 


(¥4 


(3) 


where  Z,  and  Z2  are  the  atomic  numbers  of  the  ion  and  target  atoms  and  e  is 
the  the  electronic  charge.  The  bracketed  term  is  Firsov’s  screening  function, 
with  R,  being  the  screening  radius: 


R ,  -  (9w  J/128)I/3  a0(Z'/2  +  Z\' 2 )  'V\  f4» 

where  a0  is  the  Bohr  radius  h2/me2. 

The  Firsov  potential  was  used  to  calculate  the  shadow  cone  d  s.ons  in 
this  study,  with  one  minor  correction.  It  has  been  shown  [11]  that,  if  the  above 
screening  radius  is  multiplied  by  0.8  and  used  in  the  Moliere  approxin  o 
the  Thomas-Fermi  potential,  the  ratio  of  single  to  double  scattering  for 
computer  simulation  of  2.5  keV  Ne*  scattering  from  Au  and  Ni  is  in  good 
agreement  with  experimental  results.  This  same  correction  has  been  applied  to 
the  Fnsov  potential  used  in  this  study,  so  the  potential  constant  A  is  given  by: 

A  -2.44Z,ZJ/(Z,1  :  +  Z'2'2)2  'eVA2,  (5) 

Fig.  7  shows  shadow  cones  calculated  for  9.5  keV  Ne*  on  Fe  projected  on  an 
Fe  lattice.  The  second  layer  atoms  in  this  figure  are  inside  the  shadow  cones  of 


(tiO)  ♦  -45* 


5  A 

I - 
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Fig.  7  Shadow  cones  calculaied  for  9.5  keV  Ne  ‘  scattering  from  Fe  projected  onto  an  FefllO) 
plane  for  +  -  45  °  showing  the  shadowing  of  the  second  layer  by  the  first. 
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first  layer  atoms  and  therefore  should  not  contribute  to  the  scattered  intensity. 

Using  the  shadow  cone  dimensions  given  by  eq.  (2),  the  ranges  of  <(>  and  4* 
for  which  any  atom  in  the  crystal  will  shadow  another  can  be  calculated.  This 
was  done  in  the  following  manner.  A  pair  of  atoms  was  chosen  and  the 
coordinates  of  one  with  respect  to  the  other  were  calculated  from  the  lattice 
parameter.  The  shadow  cone  was  projected  on  these  atoms  with  the  shadowing 
atom  at  the  apex  and  its  central  axis  parallel  to  the  incoming  beam  direction.  A 
plane  perpendicular  to  the  beam  direction  which  contained  the  shadowed  atom 
was  also  considered.  The  distance  /'  from  the  shadowing  atom  to  this  plane 
(see  fig.  7)  and  distance  r'  in  this  plane  from  the  shadowed  atom  to  the  central 
axis  of  the  shadow  cone  were  then  calculated.  The  distance  r'  was  then 
compared  with  r(l')  calculated  using  eq.  (2)  to  see  if  the  atom  is  inside  the 
shadow  cone.  For  each  pair  of  atoms  chosen  the  limits  of  \p  corresponding  to 
the  atom’s  being  at  the  edges  of  the  shadow  cone  were  determined  within  ±1° 
using  an  iterative  search  technique.  This  was  done  at  5°  increments  of  <J>  from 
<f»  =  0  to  <#>  =  45  ° ,  so  that  the  shadowing  and  blocking  limits  for  <f>  rotations 
were  also  determined.  These  limits  were  then  plotted  on  a  plot,  so  that 
minima  in  the  scattered  intensity  (e.g.  fig.  6)  could  be  compared  with  beam 
orientations  for  which  deeper  layers  are  shadowed. 

This  solution  gives  results  slightly  different  from  those  of  Bronckers  and  De 
Wit  [8],  who  plotted  the  shadowing  limits  as  circles  of  radius  12  = 
(180/2rr)tan~ [{r'/l')  with  r'  and  /'  determined  for  the  two  atoms  aligned 
along  the  central  axis  of  the  shadow  cone.  Bronckers  and  De  Wit  ignore  the 
variation  of  /'  when  the  shadowed  atom  is  not  on  the  central  axis  of  the 
shadow  cone;  this  is  a  good  approximation  if  /  a>  r  or  if  the  radius  of  the 
shadow  cone  is  not  a  strong  function  of  I. 

If  the  crystal  is  oriented  so  that  the  second  layer  is  exposed  to  the  beam,  it  is 
still  possible  that  the  second  layer  atoms  will  not  contribute  to  the  scattered 
intensity.  Ions  scattered  by  second  layer  atoms  can  be  deflected  or  blocked  by 
atoms  in  the  first  layer  on  the  way  out,  so  that  they  do  not  reach  the  analyzer. 
An  analysis  of  blocking  cones  has  been  done  for  experiments  in  which  the 
incident  ion  beam  was  kept  at  a  fixed  elevation  and  variations  in  intensity  were 
monitored  by  moving  the  detector  independently  [20.21].  As  shown  in  fig.  8. 
the  ion  beam  reaches  an  atom  in  the  second  (or  deeper)  layer  and  scatters  from 
it  over  a  range  of  angles,  0,.  This  atom  can  be  considered  to  be  a  point  source 
of  ions  of  various  energies  £(0,)  given  by  eq.  (1)  with  0L  =  0,.  As  the  ions  pass 
the  blocking  atom  they  will  scatter  from  it  at  different  angles.  0;.  which  will 
depend  on  £(0,)  and  their  impact  parameters  with  this  atom  which  are  also  a 
function  of  0,.  As  in  the  case  of  shadowing,  an  ion  will  scatter  through  large 
angles  if  it  approaches  the  blocking  atom  at  small  impact  parameters  and  small 
angles  for  large  impact  parameters.  The  edges  of  the  blocking  cone  correspond 
to  scattering  conditions  for  which 
d(0,  +  0,)/d0,  =0. 
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For  an  unreconstructed  Fe  surface,  ail  atoms  in  the  second  layer  are 
equivalent  so  that  shadowing  and  blocking  of  only  one  second  layer  atom  need 
to  considered.  The  regions  of  space  for  which  this  atom  is  shadowed  by 
any  of  five  nearby  first  layer  atoms  have  been  calculated,  along  with  regions 
for  which  five  equivalent  atoms  block  the  outgoing  ions.  The  atoms  chosen  are 
indicated  in  fig.  9,  which  shows  the  top  two  crystal  layers  as  seen  from  above. 
The  beam  is  assumed  to  be  incident  so  that  the  atoms  labeled  1S...5S  will 
shadow  the  atom  marked  0,  and  the  atoms  labeled  IB... SB  will  block  ions 
scattered  from  this  atom. 

The  coordinate  system  ( )  used  for  these  calculations  is  the  same  as  the 
one  used  experimentally,  and  the  range  of  <t>  shown  in  fig.  9  is  the  same  as  that 
used  for  the  intensity  plots.  Because  of  the  symmetry  of  unreconstructed  Fe 


Fe-SnlOOD  Clean 
o-1ST  Layer 
•  -  2nd  Layer 


Fig.  9.  Top  view  of  the  Fe  (001)  surface  showing  first  and  second  layer  atoms  as  light  and  dark 
circles.  Atoms  whose  shadowing  or  blocking  of  the  second  layer  atom  marked  0  have  been 
considered  are  labeled. 
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(001)  surface,  this  range  of  <f>  (45°,  from  a  (110)  azimuth  to  a  (100)  azimuth)  is 
sufficient  to  describe  all  possible  p.  Within  the  range  of  p  considered,  first 
layer  atoms  other  than  those  labeled,  and  some  second  layer  atoms,  will 
shadow  and  block  the  second  layer  atom  being  considered,  but  only  at  p  less 
than  15°  or  greater  than  75°.  These  were  chosen  as  limits  of  \p  for  which 
shadowing  and  blocking  cones  needed  to  be  calculated,  since  no  data  were 
collected  outside  this  region. 

On  a  p-p  plot  the  shadowing  limits  for  each  first  layer  atom  are  distorted 
circles,  while  the  blocking  limits  are  circles.  For  simplicity  the  shadowing  and 
blocking  limits  in  regions  where  these  circles  overlap  were  not  included.  The 
limits  of  these  combined  shadowing  and  blocking  cones  are  plotted  in  fig.  10. 
The  hatched  region  corresponds  to  p  and  p  for  which  the  second  layer  of  the 
crystal  is  shadowed  or  blocked  by  the  first.  For  comparison  the  experimental 


(H0)  (210)  (310)  OOO) 


p  (degrees) 

Fig.  10.  Plot  of  shadowing  and  blocking  limits  in  <t>  and  C  for  first  layer  Fe  atoms  on  the  2nd  layer. 
Regions  where  the  second  layer  is  hidden  are  shown  hatched.  For  comparison  an  intensity  map  on 
counts/1000)  for  Ne'  scattering  from  clean  Fe  is  also  shown. 


t 


562 


L  Marc  hut  et  a!.  /  Surface  structure  analysts  using  LEI  S  / 


clean  Fe  intensity  plot  is  also  shown  on  this  figure.  The  lowest  scattered 
intensities  are  observed  in  regions  where  the  second  layer  is  not  visible, 
according  to  the  calculation,  ft  should  be  emphasized  that  single  scattering 
only  has  been  considered  here.  Multiple  scattering  may  occur  within  the 
shaded  regions,  but  the  multiple  scattering  background  has  been  subtracted  as 
described  earlier. 

An  effect  not  satisfactorily  explained  by  shadowing  and  blocking  of  the 
second  layer  is  the  large  increase  in  scattered  intensity  from  clean  Fe  near  the 
(310)  azimuth  for  35  <  +  <  62.5.  This  might  be  explainable  in  part  as  a 
focussing  effect,  but  the  greatest  focussing  would  be  expected  near  the  edges  of 
the  shadow  cones  and  it  is  not  expected  to  persist  far  beyond  this.  Thus, 
focussing  would  not  explain  the  increased  intensity  seen  near  i/.  =  45  °,  <>  =  27 ° 
in  fig.  10.  A  similar  analysis  has  therefore  been  done  to  see  if  the  third  layer  is 


k  1 10 )  (2101  (SlO)  (1001 
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Fig.  11.  Plot  of  regions  where  third  layer  Homs  are  hidden  by  first  layer  atoms  (\\\)  and  second 
layer  atoms  (///)  as  a  function  of  ♦  and  For  comparison  an  intensity  map  for  Ne*  scattering 
from  clean  Fe  is  also  shown. 
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visible  to  the  beam  at  any  of  the  beam  orientations  studied.  As  with  the  second 
layer  analysis,  only  one  third  layer  atoms  need  be  considered,  but  shadowing 
and  blocking  of  this  atom  by  both  first  and  second  layer  atoms  must  be 
included.  The  regions  in  which  second  layer  atoms  shadow  and  block  the  third 
layer  atom  are  identical  to  those  for  which  first  layer  atoms  shadow  and  block 
the  second  layer,  as  shown  in  fig.  10.  Regions  in  which  a  third  layer  atom  is 
shadowed  and  blocked  by  first  layer  atoms  were  calculated  in  a  manner  similar 
to  that  described  above.  These  regions  are  shown  hatched  with  lines  of  positive 
slope  (\\\)  in  fig.  11.  Also  shown  are  the  shadowing  and  blocking  regions  due 
to  atoms  in  the  second  layer,  hatched  with  lines  of  negative  slope  (///).  As 
noted  previously,  these  regions  are  the  same  as  the  shadowing  and  blocking 
limits  for  first-layer  atoms  on  a  second-layer  atom.  In  fig.  11  areas  hatched 
with  positive  slope  correspond  to  regions  where  only  the  first  layer  is  visible; 
areas  hatched  only  with  negative  slope  correspond  to  regions  where  both  the 
first  and  second  layers  are  visible,  and  for  areas  with  no  hatching  the  first, 
second,  and  third  layers  should  be  visible. 

An  intensity  map  for  Ne*  scattering  from  the  clean  Fe  surface  is  also 
plotted  in  fig.  11.  The  regions  of  highest  scattered  intensity  (>  4000  counts) 
correspond  reasonably  well  to  unhatched  areas  in  this  figure  where  atoms  in 
the  first,  second,  and  third  layers  should  be  visible.  The  regions  of  lowest 
intensity  ( <  1000  counts)  coincide  with  regions  for  which  only  the  top  layer  is 
visible  (hatching  with  positive  slope). 

The  most  notable  disagreement  is  near  the  (310)  azimuth  (</>**  27°)  for 
\p  =  25°  and  ^  =  62.5°.  A  better  fit  to  the  data  would  be  observed  if  the 
interatomic  potential  were  modified  to  give  slightly  smaller  shadow  cones. 
Another  possibility  is  that  the  potential  is  correct  but  the  outermost  atomic 
planes  are  displaced  inwards  so  that  the  third  layer  is  exposed  to  the  beam  over 
a  larger  angular  range.  It  is  also  possible  that  some  multiple  scattering 
sequence  contributes  to  the  peak  at  £,/£0  —  0.47,  an  effect  not  considered  in 
this  single  scattering  analysis.  Computer  simulation  of  the  particle  trajectories 
will  be  needed  to  clarify  this  matter. 

Estimates  of  the  limits  of  error  which  could  be  introduced  by  the  neglect  of 
certain  other  effects  in  this  analysis  are  illustrated  in  fig.  12  which  shows  a 
shadow  cone  calculated  for  9.465  keV  Ne*  scattered  from  Fe  from  two  atoms 
spaced  as  in  an  Fe(110)  plane.  Firsov  estimated  that  his  potential  was  within 
±  10%  of  the  Moliire  approximation  to  the  Thomas- Fermi  potential  which  is 
generally  accepted  for  low  energies.  The  thickness  of  the  line  drawn  to 
represent  the  shadow  cone  reflects  this  uncertainty.  This  error  is  fairly  small, 
since  r{l)  varies  as  Al/i. 

Another  error  in  the  shadow  cone  analysis  arises  in  that  the  shadowed  atom 
in  fig.  12  is  assumed  to  become  visible  when  it  reaches  the  edge  of  the  shadow 
cone.  Actually,  since  the  collected  ions  scatter  through  90°  from  the  target 
atom,  they  have  approached  it  at  a  small  impact  parameter.  This  means  that 


Fig.  12.  Shadow  cone  calculated  for  9.5  keV  Ne  ’  scattering  from  Fe  showing  the  magnitude  of 
various  sources  of  error  in  this  analysts.  The  error  in  the  potential  refers  to  the  discrepancy 
between  the  dimensions  calculated  using  the  Firsov  and  Moliere-Thomas-Fermi  potentials. 


the  ions  will  scatter  from  this  atom  when  it  is  still  inside  the  shadow  cone  by 
the  present  definition.  Calculations  by  Chen  [22]  indicate  that  the  impact 
parameter  for  90°  scattering  of  9.5  keV  Ne'  on  Fe  is  about  0.1  A.  This 
distance  is  shown  as  an  arrow  at  the  shadowed  atom  in  fig.  12,  and  the  result  is 
that  the  shadow  cone  calculated  is  2°  or  3°  too  large. 

There  will  be  some  uncertainty  about  the  positions  of  both  the  shadowing 
and  shadowed  atom  due  to  vibrations.  The  error  bars  in  fig.  12  indicate 
expected  rms  vibration  amplitudes  of  the  atoms.  These  were  taken  from 
calculations  of  Gupta  and  Hemkar  [23]  which  were  based  on  the  variation  of 
the  Debye- Waller  factor  with  temperature  [24],  and  theoretical  predictions  of 
Dobrzynski  and  Masri  [25]  from  the  elastic  constants  of  Fe.  The  effect  of 
vibrations  will  be  to  superimpose  a  Gaussian  distribution  on  the  step  function 
that  might  otherwise  to  be  expected  to  describe  the  change  in  intensity  near  the 
edge  of  a  shadow  cone. 

Another  change  in  scattered  intensity  near  the  edge  of  the  shadow  cone  will 
be  enhancement  of  the  intensity  just  outside  the  shadow  cone  due  to  focussing. 
This  occurs  as  many  incoming  ions  are  deflected  through  small  angles  away 
from  the  shadowed  region  and  into  the  nearby  region  just  outside  the  shadow 
cone,  so  that  an  atom  in  this  region  receives  a  higher  than  average  flux.  The 
expected  width  of  the  region  of  enhanced  intensity  is  indicated  by  the  arrows 
just  outside  the  shadow  cone  in  fig.  12  and  is  taken  from  De  Wit  and 
Bronckers  [7],  who  calculated  the  distribution  of  5  keV  Ne*  ions  scattered  by  a 
Cu  atom  using  an  A/r2  potential  similar  to  the  one  used  here. 

This  analysis  of  the  clean  Fe(001 )  surface  demonstrates  that  LEIS  can  be 
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used  for  layer  by  layer  chemical  analysis  of  unreconstructed  surfaces  since 
sample  orientations  at  which  different  layers  are  exposed  to  the  ion  beam  can 
easily  be  predicted.  The  focussing  factors  for  each  individual  layer  should  be 
determined  if  possible  when  performing  such  an  analysis.  A  similar  analysis 
[26]  has  been  done  for  an  ordered  Cu3Au  (001)  surface  and  good  agreement 
between  measured  and  predicted  intensity  variations  with  sample  orientation 
were  also  observed.  Another  possible  application  of  the  shadow  cone  analysis 
is  to  test  various  models  of  reconstructed  surfaces. 
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